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INTRODUCTION 

The equation of state (EOS) of superdense matter in 
neutron star cores is still a mystery. It is not clear if it 
is soft, moderate or stiff; if the matter contains nucle- 
ons/hyperons, or exotic components. In the absence of 
good practical theory of supranuclear matter the prob- 
lem cannot be solved on purely theoretical basis, but it 
can be solved by comparing theoretical models with ob- 
servations of neutron stars. The attempts to solve this 
long-standing problem by different methods are numer- 
ous (e.g., Refs. [1, 2]). Here we discuss current results 
obtained from studies of cooling isolated neutron stars. 

The first papers on neutron star cooling appeared with 
an advent of X-ray astronomy, before the discovery of 
neutron stars. Their authors tried to prove that neutron 
stars cool not too fast and can be discovered as sources 
of thermal surface X-ray radiation. The first estimates of 
thermal emission from cooling neutron stars were most 
probably done by Stabler [3] in 1960. Four years later 
Chiu [4] made similar estimates and analyzed the possi- 
bility to discover neutron stars from their thermal emis- 
sion. First, simplified calculations of neutron star cooling 
were done in 1964 and 1965 [5, 6, 7]. The foundation of 
the strict cooling theory was laid in 1966 by Tsuruta and 
Cameron [8], one year before the discovery of pulsars. 
We review the current state of the cooling theory. More 
details can be found in recent review papers [9, 10]. 

COOLING NEUTRON STARS 

Neutron stars are thought to consist of a thin crust and 
a massive core (e.g., Ref. [1]). The crust constitutes a 



few per cent of the star's mass and is about ^ 1 km 
thick. The mass density p at the crust-core boundary is 
« Po/2, where po « 2.8 x 10^"* g cm^^ is the density of 
saturated nuclear matter. The crust is further divided into 
the outer and inner crust. The outer crust (more generally, 
the outer envelope; where p is below the neutron drip 
density < 4 x lO'^ g cm^^) consists of atomic nuclei 
and strongly degenerate electrons. In the inner crust, free 
neutrons appear in addition to the electrons and nuclei; 
these neutrons can be in superfluid state. 

The neutron star core can be divided into the outer and 
inner core. The outer core extends to p ^ (2 — 3)po and 
consists of neutrons (n), with an admixture of protons 
(p), electrons (e) and possibly muons. All these particles 
are strongly degenerate. The inner core extends to the 
stellar center [to ~ (10 — 20)po in most massive stars]. 
Its composition is very uncertain. It may be the same 
as in the outer core or essentially different. In particular, 
hyperons may appear there in addition to nucleons. An- 
other possibility is the appearance of exotic matter (pion 
or kaon condensates or strange quark matter or mixed 
phases, as reviewed, e.g., in [1]). Nucleons, hyperons and 
quarks can be in superfluid state. Physical properties of 
matter at p < po are more or less restricted by nuclear 
physics data, but at higher p they are uncertain. Typi- 
cal baryon chemical potentials ^ 500 MeV in a neutron 
star core are most difficult for rigorous microscopic cal- 
culations. Note a lively discussion in the literature that 
compact stars (or some of them) can be not neutron stars 
but strange stars (e.g., Refs. [1, 2, 11]) built entirely or 
almost entirely of strange quark matter. 

Neutron stars are born in supernova explosions with 
high internal temperature T ^ lO" K, but gradually cool 
down. In ^30 s after the birth a star becomes trans- 
parent for neutrinos generated in its interiors and trans- 
forms from a proto-neutron star (e.g., [12]) to an ordinary 
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TABLE 1. Main slow neutrino emission in nucleon matter 



Process 






Qs [erg cm ^ s ^] 


Modified Urea 


n + N- 


-^p + N + e+v p + N + e-^n+N+v 


lO^" -3x10^1 


Bremsstrahlung 


N + N 


^N + N+v + v 


IQiy - io20 



neutron star whose EOS is almost temperature indepen- 
dent (except near the very surface). At the later neutrino- 
transparent stage the star cools via neutrino emission 
from its interior and via heat transport to the surface and 
subsequent photon thermal emission. 

Neutron star cooling is described by general relativis- 
tic equations [13] of heat diffusion inside the star with 
neutrino energy sources and surface photon emission. 
The solution gives the distribution of the temperature T 
inside the star versus time t, and the effective surface 
temperature Tfi(t). To simplify computations, one divides 
[14] the neutron star into the interior region and the out- 
ermost heat-blanketing envelope [extended to densities 
P < Ph ^ (10"^ - 10") g cm"^]. The thermal structure 
of the blanketing envelope is studied separately in the 
stationary, plane-parallel approximation, that relates 
to the temperature Tj, at p = Pb- The diffusion equations 
are then solved in the interior (p > pb). 

The thermal photon luminosity of the star is Ly = 
AnaR^T^{t), where R is the circumferential stellar ra- 
dius. Both, Ly and T^, refer to a locally-flat reference 
frame on the stellar surface. A distant observer detects 
the "apparent" luminosity L" = Ly{l — r^/R) and the 
"apparent" effective temperature = y'l — r^/R, 
where = IGM/c^ is the Schwarzschild radius and M is 
the gravitational neutron star mass. If the surface temper- 
ature distribution is anisotropic (for instance, owing to a 
strong magnetic field), Ly is determined by a properly 
averaged surface temperature (e.g., Ref. [15]). 

The theory gives cooling curves, T^{t), and predicts 
three main cooling stages. The initial stage of internal 
thermal relaxation lasts f < 10 — 50 yr. At this stage, the 
neutron star crust stays hotter than the core and thermally 
decoupled from it (because of much stronger neutrino 
cooling in the core). The surface temperature is then in- 
sensitive to the physics of the core but strongly depends 
on physical properties of the crust [16, 17]. Since no neu- 
tron star has been observed at this stage, we wiU not dis- 
cuss it in detail. The next stage of neutrino cooling with 
isothermal interior lasts f < 10^ — 10^ yr. At this stage the 
neutrino luminosity Ly ^ Ly\ thermal conduction is high 
and makes the stellar interior isothermal, with the main 
temperature gradient located in the heat blanketing en- 
velope. The cooling is mostly regulated by a strong neu- 
trino emission from the core. The surface temperature re- 
sponds to the core cooling and depends on properties of 
superdense core. All cooling isolated neutron stars whose 



thermal surface radiation has been detected seem to be at 
this neutrino cooling stage (or its very end). Therefore, 
we focus our attention on this stage. At the next photon 
cooling stage the star is cold; its neutrino emission dies 
out (Ly <C Ly) and the cooling is governed by photon sur- 
face emission. 

Soon after the coohng starts (in minutes to a year, 
depending on the internal structure), the core temperature 
drops to r ^ 10^ K. For this T, the internal thermal 
energy of the star is ~ lO'*^ erg. The cooling theory 
shows how this heat emerges from the star. 

NEUTRINO ElMISSION IMECHANISMS 

Let us summarize the main neutrino emission mecha- 
nisms in the neutron star core. They are strongly affected 
by baryon superfluidity. More details can be found, e.g., 
inRefs. [18, 19, 20]. 

Neutrino emission in nonsuperfluid cores. The ma- 
jor neutrino mechanisms in nucleon matter of the outer 
core (p < 2 po) are modified Urea process and nucleon- 
nucleon bremsstrahlung. They are listed in Table 1 (from 
Ref. [9]), where N denotes a nucleon (n or p). They 
are relatively weak and produce slow neutrino cooling. 
The modified Urea process differs from its direct Urea 
progenitor, described below, by an additional nucleon- 
spectator that is required to satisfy momentum conserva- 
tion of reacting particles. All reactions involving elec- 
trons can involve muons instead (if available in dense 
matter). 

At higher p, in the inner core, neutrino emission can 
be strongly enhanced by new mechanisms (Table 2, from 
[9]). The enhancement level greatly depends on the EOS 
and composition of superdense matter that is most im- 
portant for the cooling problem. The strongest enhance- 
ment is provided by direct Urea process [21, 22] in nu- 
cleon (or nucleon-hyperon) matter. It is a sequence of 
two reactions (a beta decay and beta capture) forbidden 
in the outer core by momentum conservation. It is al- 
lowed in the inner core for the matter with rather high 
( > 11-13%) proton fraction (for EOSs with large sym- 
metry energy of nuclear matter). If it is forbidden but the 
matter contains pion condensate, neutrino emission is en- 
hanced (although weaker) by direct-Urca-type reactions 
involving quasinucleons N (superpositions of n and p) 
in pion-condensed matter. If pion condensate is absent. 
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TABLE 2. Leading processes of fast neutrino emission in nucleon matter and 
three models of exotic matter 



Model Process Qf [erg cm ^ s ^ ] 



Nucleon matter 


n 


^p + e + 


V 


p + e ^ 


n + V 


1q26_ 


3 X 102^ 


Pion condensate 


N 


-^N + e + 


V 


N + e^ 


► N + V 


1023 


_1026 


Kaon condensate 


B 


-^B + e + 


V 


B + e^ 


B + v 


1023 


-1024 


Quark matter 


d 


^u + e + 


V 


u + e —> 


d+ V 


1023 


-1024 



but kaon condensate available, neutrino emission is en- 
hanced (even more weaker) by direct-Urca-Uke process 
involving baryonic quasiparticles B in kaon condensed 
matter. Nearly the same enhancement is expected due to 
the direct Urea process involving d and u quarks in quark 
matter. The processes in the inner core (Table 2) can am- 
plify the neutrino emission by 2-7 orders of magnitude 
and lead to fast cooling. It is also possible that the neu- 
trino emission in the inner core is not enhanced. 

The emissivity Qv{p,T) of slow and fast neutrino 
processes in nonsuperfluid matter can be written as 

Gsiow = QsTg, Gfast = 2f 5 ( 1 ) 

where = r/(10^ K), while Qs and Qf are slowly vary- 
ing functions of p (Tables 1 and 2). Thus, the neutrino 
luminosities of massive and low-mass stars (with and 
without iimer core) can be very different. For instance, 
in a young massive star, where the direct Urea is open 
and the core temperature is T = 10^ K, Ly is as huge 
as ^ 10^^ erg s^'. In a low-mass star at the same T, Ly 
would be ~ 7 orders of magnitude lower. In both cases Ly 
rapidly decreases when the star cools. For instance, one 
has Ly ~ 10^'* erg s~^ for a low-mass star at ? ~ 10 kyr. 

Neutrino emission in superfluid cores. Nucleons, hy- 
perons, and quarks in dense matter can be in superfluid 
state (e.g., [23], also see [9, 19] for references). This su- 
perfluidity occurs via Cooper pairing of particles owing 
to an attractive component of their interaction (with the 
appearance of a gap in the particle energy spectrum near 
the Fermi level). Systematic simulations of cooling su- 
perfluid neutron stars were triggered by a remarkable pa- 
per of Page and Applegate [24]. Superfluidity of any par- 
ticles is characterized by its own density dependent criti- 
cal temperature rc(p). Microscopic calculations of rc(p) 
are extremely model dependent and give a large scatter 
of critical temperatures. Let us mention several general 
features. 

Neutron pairing in the spin singlet state with zero an- 
gular momentum (^Sq) occurs in the inner neutron star 
crust (for free neutrons) and dies out in the core be- 
cause singlet-state nuclear attraction turns into repulsion 
near the crust-core interface. Neutron pairing in the core 
may occur in a triplet state with unit angular momen- 
tum, (coupled to a spin triplet state with three units 



of angular momentum, 3F2). Cooper pairing of other 
particles in the core can occur either in a singlet state 
or in a triple state. Baryon superfluidity is affected by 
the presence of pion or kaon condensate (e.g., [25] and 
references therein). Calculated values Tc{p) range from 
~ 10* K to a fewxlO'*^ K and vanish at supranuclear 
densities (where attractive pairing interaction becomes 
inefficient). 

A special case is presented by color superconductivity 
[26] owing to very strong pairing of unlike quarks in 
quark matter, where Tc (p ) can be as high as ^ 5 x 1 " K. 

Any baryon superfluidity reduces neutrino processes 
involving these baryons (Tables 1 and 2) due to a gap in 
the baryon energy spectrum. At T ^ Tc, the reduction 
is exponentially strong. For instance, the powerful direct 
Urea process can be formally open in the inner core, 
but completely suppressed by a strong superfluidity of 
neutrons of protons. 

In addition, when T falls below T^, superfluidity initi- 
ates a specific neutrino emission due to Cooper pairing 
of baryons [27] which enhances neutrino cooling. This 
mechanism can enhance the neutrino luminosity of the 
star by a factor of < 30 — 100 over the modified Urea 
level [28, 29]. 

Superfluidity affects also neutron star heat capacity 
[19] and thermal conductivity [30, 31] but these effects 
are less strong. 

OBSERVATIONS 

We will compare cooling theory with observations of sur- 
face thermal radiation of isolated neutron stars. The cur- 
rent state of the observations is discussed, e.g., in Refs. 
[32, 33, 34]. In Table 3 we list 16 isolated neutron stars 
whose effective surface temperatures have been mea- 
sured or constrained. For brevity, the stars are numbered, 
and these numbers are used in the text and figures. 

The data include two neutron stars (1 and 2, the pul- 
sars Crab and J0205+6449) in historical supernova rem- 
nants (SNRs); the famous Vela pulsar (7) and the similar 
PSR B 1706-44 (8); the "pulsar-twins" J0538+2817 and 
B2334H-61 (9 and 10); the "three musketeers" [Geminga, 
PSR B1055-52, and PSR B0656+14 (12, 14, and 11)]; 
one compact central object in SNR (RX J0822-4300, 
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TABLE 3. Observational limits on surface temperatures of isolated neutron stars 



Number 


Source 


/ [kyr] 


h LMKJ 




Model 


Kei. 


1 


PSRB0531+21 (Crab) 


1 


<2.0 


99.8% 


BB 


[37] 


2 


PSR J0205+6449 (in 3C 58) 


0.82-5.4 


<1.02 


99.8% 


BB 


[38] 


3 


PSRJl 119-6127 


~ 1.6 


f» 1.2 


- 


mHA 


[32] 


4 


RX J0822-4300 (in Pup A) 


2-5 


1.6-1.9 


90% 


HA 


[39] 


5 


PSRJ1357-6429 


~7.3 


fa 0.766 


- 


mHA 


[40] 


6 


RX J0007. 0+7303 (in CTA 1) 


10-30 


<0.66 


- 


BB 


[41] 


7 


PSR B0833-45 (Vela) 


11-25 


0.68 ±0.03 


68% 


mHA 


[42] 


8 


PSR B 1706-44 


~17 


0.82+°;»i 


68% 


mHA 


[43] 


9 


PSR J0538+2817 


30±4 


~ 0.87 




mHA 


[44] 


10 


PSRB2334+61 


~41 


~0.69 




mHA 


[32] 


11 


PSRB0656+14 


~110 


0.91±0.05 


90% 


BB 


[45] 


12 


PSRB0633+1748 (Geminga) 


~340 


~0.5 




BB 


[46] 


13 


RX J 18564-3754 


~500 


0.434 ±0.003 


68% 


mHA* 


[47] 


14 


PSR B 1055-52 


~540 


~0.75 




BB 


[48] 


15 


PSR J2043+2740 


~ 1200 


~0.44 




mHA 


[32] 


16 


RX J07204-3125 


~ 1300 


~0.51 




HA* 


[49] 



source 4); the compact source RX J0007.0+7303 at the 
center of the SNR CTA 1 (6); two "dim" ("truly iso- 
lated") neutron stars RX J 1856.4-3754 and RX J0720.4- 
3125 (13 and 16); two young and energetic pulsars 
Jl 1 19-6127 and J1357-6429 (3 and 5); and the old pul- 
sar J2043+2740 (15). 

The ages and effective surface temperatures of many 
sources are rather uncertain. For the sources mentioned 
in [29, 35, 36] the choice of f and is mainly the 
same as in these references. The ages of other sources 
are pulsar spindown ages, and the errobars of t and 
are chosen in the same way as in [29, 35, 36]. We have 
enlarged a possible age range of the source 2 to 5.4 kyr, 
the pulsar spindown age, because this source can be acci- 
dentally projected onto the SNR 3C 58 (Yu.A. Shibanov, 
private communication, 2007). The values of r^" are in- 
ferred from observations (reported in the indicated ref- 
erences) assuming either black-body (BB) spectrum or a 
hydroden atmosphere model (magnetic or nonmagnetic 
one, mHA or HA). For the sources 13 and 16, is de- 
termined [47, 49] using models of hydrogen atmospheres 
of finite depth (mHA* and HA*, respectively). One im- 
portant source, IE 1207.4-5209, is not included in the 
table because of the problems to interpret its spectrum 
(although modeling of such spectra is progressing [50]). 

BASIC COOLING CURVES 

For illustration, we will mainly use models of neutron 
stars whose cores contain neutrons, protons and elec- 
trons and have a stiff phenomenological EOS proposed 
in Ref. [51] (model 1 for the symmetry energy and the 
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logp (g cm-=') log p (g cm-=') 

FIGURE 1. (a) Superfluid transition temperature versus 
density for two models (Ip and 2p) of proton superfluidity in 
the neutron star core, (b) Neutrino emissivity profiles in the 
core at r = 3 X 10** K for nonsuperfluid matter (noSF) and for 
matter with superfluid protons (models Ip or 2p). 



bulk energy model that gives the compression modulus 
of saturated nuclear matter K = 240 IVleV). The most 
massive stable neutron star, for this EOS, has the grav- 
itational mass Mmax = 1.977M0 and the central density 
Pc = 2.578 X lO'^ g cm^^; the direct Urea process is al- 
lowed at p > Pd = 7.851 X lO''* g cm^-' (in stars with 
M>Md = I.358M0). 

Neutrino emissivity profiles in a stellar core at T = 
3 X 10^ K are shown in Fig. lb. The vertical dotted 
line shows the direct Urea threshold, p^. The thick line 
(NoSF) is for nonsuperfluid matter At p < po the neu- 
trino emission is determined by the modified Urea pro- 
cess. When p increases above po, the emissivity jumps 
by 7 orders of magnitude due to the direct Urea onset. 

Fig. 2a presents the cooling curve for a 1.3M0 non- 
superfluid star compared with the data. The curve is a 
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log t (yrs) log t (yrs) log t (yrs) log t (yrs) 

FIGURE 2. Observational limits on surface temperatures of isolated neutron stars compared to theoretical cooling curves. 

(a) Cooling curve for a nonsuperfluid 13Mq, star. Any other panel shows 88 curves from M = 1.1 to 1.97M0 with step O.OIM©: 

(b) No superfluidity, sharp direct Urea threshold; (c) proton superfluidity Ip and (d) 2p which broaden the threshold. 
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FIGURE 3. Very schematic view of (a) neutrino emissivity 
profile in a neutron star core and (b) neutrino luminosity versus 
star's mass at a fixed core temperature T. Units are arbitrary. 



typical example of slow cooling via the modified Urea 
process. It is actually universal, being almost the same 
for all stars with 1 . 1 Mq <M< M-q, and for a wide class 
of EOSs [24]. This universal curve goes through the scat- 
ter of observational points but cannot explain the sources 
which are hottest and coolest for their age. These sources 
seem to have neutrino luminosities lower and higher than 
the modified Urea level. 

Now consider a set of nonsuperfluid neutron star mod- 
els of masses from 1.1 to \ .91 Mq, with the mass step of 
0.01 Mq (Fig. 2b). The first 26 models (M < 1.35M0 < 
Md) belong to the class of low-mass stars and show the 
universal slow cooling behavior. The next 1.36M0 star 
has a small central kernel with the enhanced neutrino 
emission and shows much faster cooling (is much colder 
than all observed sources). All other more massive stars 
cool even faster via the direct Urea process and belong 
to the family of cold massive rapidly cooling stars. Evi- 
dently, these cooling curves cannot explain the data. 



BASIC CONCEPT 

A majority of realistic cooling scenarios, constructed up 
to now, satisfy the following basic phenomenological 
concept: 

(1) Low-mass stars undergo slow neutrino cooling. 

(2) Massive stars cool much faster via neutrino emis- 
sion from their inner cores. 

(3) There is a family of medium-mass stars whose 
cooling is intermediate between the slow and fast ones. 

Accordingly, the density profile of the neutrino emis- 
sivity in a neutron star core can look like that plotted in 
Fig. 3a; pi and p2 mark the density range, where the slow 
cooling transforms into the fast one. This 2v(p) profile 
translates into a characteristic dependence of the neutrino 
luminosity Ly on the stellar mass shown in Fig. 3b; the 
transition between the slow cooling to the fast one takes 
place in the mass range from M\ to M2. 

Thus, at the present stage of investigation, the cooling 
theory has potential to test four main parameters of the 
neutrino emissivity function 2v(p) or the luminosity 
function Ly {M) . They specify the lower and upper levels 
of neutrino emission (either Qsiow and Qfast or Lsiow and 
Lfast) and the position of the transition zone between the 
slow and fast emission (either pi and p2 or Mi and M2). 



EXAMPLES OF PHYSICAL MODELS 

Two examples [52] of successful realization of the above 
scheme are given in Figs. 2c and d. The figures show the 
cooling curves for neutron stars ofM = (1.1 — 1.97) Mq 
with strong proton superfluidity (Ip or 2p) in the core 
(but with normal neutrons). The density profiles rc(p) of 
the proton critical temperature are shown in Fig. la; they 
are phenomenological and can be regarded as illustra- 
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tive. In the outer core, proton superfluidity is strong and 
suppresses the modified Urea process. Accordingly, the 
slow neutrino emission of low-mass stars is determined 
by weaker neutrino bremsstrahlung in neutron-neutron 
collisions. This rises the surface temperature of low-mass 
stars and allows one to explain observations of hottest 
sources. 

In addition, proton superfluidity penetrates into the in- 
ner core (p > Pd) and suppresses the direct Urea pro- 
cess within the penetration depth. It broadens the di- 
rect Urea threshold and realizes a smooth transition from 
slowly to rapidly cooling neutron stars with the growth 
of M. Superfluidity 2p penetrates deeper than Ip which 
widens the range of masses of medium-mass stars. One 
can now explain observations of all the sources, includ- 
ing the coldest ones. Moreover, these scenarios predict 
the existence of very cold massive neutron stars (which 
have not been observed so far). 

If the EOS and the critical temperature Tc{p) in the 
stellar core were known, one would be able to "weigh" 
neutron stars [52]. For instance, the mass of the Vela pul- 
sar would be 1.47M0 for superfluidity Ip and I.6IM0 
for 2p. Unfortunately, neither EOS nor Tc{p) are known 
and the "weighing" procedure is ambiguous. Neverthe- 
less, as a rule, the ambiguity does not destroy the mass 
ordering of cooling stars. For instance, the Vela pulsar 
is expected to be more massive than RX J0822-4300. 
Moreover, one can invert the scheme and state that had 
the mass of one or several cooling neutron stars been 
measured, one could constrain the EOS and Tc{p). 

One can construct many other models of cooling neu- 
tron stars consistent with the data. One can take other 
EOSs in the core [with other direct Urea thresholds po 
and functions Tc{p)]. For instance, one can assume the 
presence of strong neutron superfluidity and normal pro- 
tons (because cooling curves are nearly symmetric [53] 
with respect to exchanging neutron and proton Tc{p) pro- 
files). 

The direct Urea threshold can also be broadened by 
nuclear effects (e.g., of pion polarization; see [54] and 
references therein) or by superstrong magnetic fields 
> 3 X 10'^ G [55] in a nonsuperfluid neutron star core. 

LEVELS OF SLOW AND FAST 
NEUTRINO EMISSION 

After several examples we can turn to a general analysis. 

One can show that it is sufficient to lower the neu- 
trino luminosity due to the modified Urea process in a 
low-mass star by a factor of ~ 30 — 100 to explain ob- 
servations of stars hottest for their age. This low neu- 
trino luminosity can be produced by nucleon-nucleon 
bremsstrahlung. The modified Urea process can be re- 



duced by nucleon superfluidity, as in the above examples 
(Figs. 2c and d). 

The level of fast neutrino emission in high-mass stars 
is less certain. To analyze it, let us consider (Fig. 4) 
three typical levels appropriate to the direct Urea process 
in nucleon/hyperon matter; to pion-condensed matter; 
and kaon-condensed matter. A schematic picture of the 
neutrino emissivity profiles through the steUar core for 
these scenarios is presented in Fig. 4a. 

Fig. 4b shows the ranges of (hatched regions) 
which can be explained by the theory assuming the sce- 
nario with the direct Urea process. The upper curve 
shows cooling of low-mass stars via neutron-neutron 
bremsstrahlung neutrino emission (as in Figs. 2c and d). 
The lower curve presents fast cooling of a maximum- 
mass star with direct Urea process open in the inner core. 
Any value of r~ between these two curves can be ex- 
plained by a cooling of a star with some intermediate 
mass. As discussed above, this scenario explains all the 
data and predicts the existence of very cold neutron stars. 

Fig. 4c refers to the scenario with pion-condensation. 
The upper cooling curve is the same as in Fig. 4b (be- 
cause low-mass stars have no inner cores), but the lower 
curve goes essentially higher (because the neutrino emis- 
sion due to pion condensation is lower than due to direct 
Urea process). This scenario also explains all the sources 
and predicts the existence of cold stars, but not so cold as 
in the direct Urea case. 

Fig. 4d is a similar plot for the scenario with kaon- 
condensation in the inner core. High-mass stars cool 
slower than in the previous scenarios, but all the data are 
still explained and a population of stars slightly colder 
than the observed ones is predicted. 

Thus, all three scenarios are currently consistent with 
the observations, but predict different families of cold 
stars. Similar conclusions have been made by a number 
of authors (see, e.g., [9, 10] and references therein). Cal- 
culations show that the the neutrino luminosity in the in- 
ner core should be w30-100 times higher than the modi- 
fied Urea luminosity to explain all the sources and leave 
no space for neutron stars colder than the observed ones. 
A discovery of very cold neutron stars would be crucial 
to firmly constrain the enhanced neutrino emission level. 

Unfortunately, current observations give no reliable 
constraints on the position of the transition layer between 
slowly and fastly neutrino emitting regions in the core 
(on pi and p2 or Mi and M2 in Fig. 3). 

OTHER ASPECTS 

We do not discuss in detail many aspects of the cooling 
theory but mention them here. 

First, cooling can also be regulated by heat conduction 
in the neutron star crust. The thermal conductivity in the 
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FIGURE 4. (a) Schematic density dependence of tiie neutrino emissivity in a neutron star core at T = 10^ K assuming very slow 
neutrino emission in the outer core and three scenarios of fast emission in the inner core, (b), (c), and (d): Ranges of (hatched) 
for the three types of fast emission compared with observations (see text for details). 



crust can be affected by composition of the matter (e.g., 
by the presence of light elements in accreted surface 
layers) and magnetic fields (in the outer and inner crust; 
see, e.g., [9, 10, 56] and references therein). However, 
these cooling regulators are usually not as strong as 
neutrino emission from the core. 

The effects of superfluidity can be more sophisticated 
than described above. Taking different EOSs [different 
direct Urea thresholds and Tc (p ) profiles] one can obtain 
a complicated zoo of families of cooling neutron stars 
[35, 57]. Neutrino emission due to Cooper pairing of 
baryons can affect the cooling in many ways as reviewed 
in Refs. [9, 10]. In particular, a moderately strong neu- 
tron superfluidity in the nucleon core can enhance the 
neutrino luminosity over the modified Urea level by a 
factor of ~30, allowing one to explain the coldest ob- 
served sources by neutrino emission from the nucleon 
(non-exotic) core without invoking the direct Urea pro- 
cess. This leads to the scenario of "minimal cooling" 
[28, 29]. 

Cooling of neutron stars can be accompanied by some 
reheating associated, for instance, with the dissipation of 
differential rotation, deviation from beta-equilibrium or 
decay of magnetic fields in neutron star interiors. The 
magnetic reheating can be especially strong in magnetars 
(see below). Otherwise these reheating mechanisms are 
expected to be significant in old neutron stars (at the 
photon cooling stage). They are reviewed in [9, 10]. 

Finally, we have not discussed cooling of strange stars 
and neutron stars with quark cores. A comprehensive 
review is given in [10]. 

CONNECTIONS 

With the spectacular progress of neutron star observa- 
tions, cooling problem becomes more connected to other 



aspects of neutron star physics. 

First, cooling of isolated neutron stars is closely re- 
lated to deep crustal heating of transiently accreting neu- 
tron stars in soft X-ray transients (SXTs). The matter ac- 
creted on a neutron star becomes eventually compressed 
by the weight of newly accreted material and sinks thus 
into the crust. Sinking into the deep crust is accompanied 
by nuclear transformations [58, 59, 60, 61] and associ- 
ated heating of the entire star (1.5 — 2 MeV per one ac- 
creted nucleon). This deep crustal heating can keep ther- 
mally inertial stars warm [62] even in quiescent states of 
SXTs when accretion is switched off. It is likely respon- 
sible for thermal emission observed from some SXTs in 
quiescent states (e.g., Refs. [63, 64]). 

One can distinguish two heating regimes. As a rule, 
accretion episodes are not too long (weeks-months) and 
do not violate isothermality of internal layers. These stars 
reach thermal quasiequilibrium determined by the mass 
accretion rates (M) averaged over their global thermal 
relaxation time 100 — 1000 yr). Observations of these 
objects test essentially the same physics of neutron star 
interior as observations of cooling isolated neutron stars 
(e.g., [65, 66, 67]). The analysis of the data on SXTs 
[63, 64, 65] generally agrees with the above analysis 
of cooling neutron stars, but with important exceptions. 
The upper limits on in quiescent states of the tran- 
siently accreting millisecond pulsar SAX J1808. 4-3658 
and possibly the SXT IH 1905H-000 are very low, giving 
an example of very cold neutron stars (missing so far in 
the data on cooling neutron stars). If true, these neutron 
stars undergo strong neutrino cooling via direct Urea pro- 
cess. However, in view of uncertainties associated with 
the deep crustal heating hypothesis, determination of Ts 
and (M), these results should be taken with caution. 

Some SXTs (KS 1731-260 and MXB 1659-29, e.g., 
[68]) show long accretion episodes (a few years and 
more) in which the crust may become much warmer than 
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the core. After accretion stops, the crust cools and ther- 
mally equilibrates with the core. This crust-core relax- 
ation is reflected in the relaxation of the surface temper- 
ature Ts{t), that is observed and can give important infor- 
mation on the crust and core physics [68, 69, 70, 71]. 

Coohng theory is also employed to study thermal 
states of magnetars. It seems one needs to introduce a 
reheating (most probably produced by strong magne- 
tar magnetic fields) to explain high observed X-ray lu- 
minosities of magnetars. One can study quasistationary 
thermal states of magnetars (e..g., [72]) and afterburst re- 
laxation (e.g., [73]). 

In addition, cooling theory can be used to analyze the 
propagation of thermal waves in neutron stars during X- 
ray bursts and superbursts (e.g., [74, 75|) and to explore 
coohng of vibrating neutron stars [76]. 

Since the data are becoming numerous it is possible to 
combine cooling theory with neutron star statistics [77]. 

CONCLUSIONS 

The theory of cooling neutron stars of ages 10^ — 10'^ 
yr mostly tests the neutrino emission properties of the 
neutron star core. Its main results are as follows. 

(1) Neutrino emission in the outer core (i.e., in the core 
of a low-mass star) is a factor of 30-100 lower than the 
modified Urea emission in a nonsuperfluid star. 

(2) Neutrino emission in the inner core (of a massive 
star) is at least a factor of 30-100 higher than the mod- 
ified Urea emission. It can be enhanced by direct Urea 
process in nucleon/hyperon inner core or by the presence 
of pion or kaon condensate, or quark matter. 

(3) The scenario with open direct Urea process pre- 
dicts the existence (Fig. 4) of massive isolated neutron 
stars which are much colder than those observed now. 
In the scenario with pion condensate, the massive stars 
should be warmer (than those with open direct Urea) but 
colder than the observed ones. In the scenario with kaon 
condensate the massive stars should be even warmer but 
slightly colder than the observed sources. A discovery of 
cold cooling neutron stars would be crucial to constrain 
the level of enhanced neutrino emission in the inner core. 

(4) Observations of cooling neutron stars can be an- 
alyzed together with observations of SXTs in quiescent 
states. The data on SXTs indicate the existence of very 
cold neutron stars (first of all, SAX J1808.4-3658) which 
cool via direct Urea process, but the data and interpreta- 
tion require additional confirmation. 

(5) A transition from slow neutrino emission in the 
outer core to enhanced emission in the inner core has 
to be smooth. Current observations of cooUng neutron 
stars and SXTs do not constrain the parameters of this 
transition. A firm measurement of masses of cooling or 
accreting stars would help to impose such constraints. 



(6) New observations and reliable practical theories 
of dense matter are vitally important to tune the cooling 
theory as an instrument for exploring physical properties 
of neutron star interiors and neutron star parameters. The 
tuning will imply a careful analysis of many cooUng 
regulators. 
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